Draft version August 19, 2010 

Preprint typeset using IAT^X style cmulatcapj v. 03/07/07 



THE [O II] A3727 LUMINOSITY FUNCTION AT Z ~ 1 
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ABSTRACT 

We measure the evolution of the [O II] A3727 luminosity function at 0.75 < z < 1.45 using high- 
resolution spectroscopy of ~ 14, 000 galaxies observed by the DEEP2 galaxy redshift survey. We 
find that brighter than L[on] = 10 42 erg s _1 the luminosity function is well- represented by a power 
law dN/dL cx L a with slope a ~ —3. The number density of [O II]-emitting galaxies above this 
luminosity declines by a factor of > 2.5 between z ~ 1.35 and z ~ 0.84. In the limit of no number- 
density evolution, the characteristic [O II] luminosity, £jpii]> defined as the luminosity where the space 

density equals 10~ 35 dex -1 Mpc~ 3 , declines by a factor of ~ 1.8 over the same redshift interval. 
Assuming that £[oii] is proportional to the star-formation rate (SFR), and negligible change in the 
typical dust attenuation in galaxies at fixed [O II] luminosity, the measured decline in £f ii] implies 
a ~ 25% per Gyr decrease in the amount of star formation in galaxies during this epoch. Adopting 
a faint-end power-law slope of —1.3 ± 0.2, we derive the comoving SFR density in four redshift bins 
centered around z ~ 1 by integrating the observed [O II] luminosity function using a local, empirical 
calibration between £[oii] and SFR, which statistically accounts for variations in dust attenuation 
and mctallicity among galaxies. We find that our estimate of the SFR density at z ~ 1 is consistent 
with previous measurements based on a variety of independent SFR indicators. 
Subject headings: galaxies: evolution — galaxies: luminosity function — stars: formation 



1. INTRODUCTION 

Measuring the comoving space density of the star for- 
mation rate (SFR) as a function of cosmic epoch is one 
of the key issues concerning the study of galaxy forma- 
tion and evolution. The current picture is that the star 
formation peaked at z ~ 1 — 3, and then declined by 
roug hly an order- of-ma gnitude to the present day (see , 
e.g.. iMadau et al.lll99l iLillv et al J 11999: lHopkinsil2004l: 
iHopkins fc Beacomll2006t) . 

Among the most direct indicators of the instanta- 
neous i_SFl^Jiigal_a2dBS is_the Ha A6563 recombination 
line (|Kennicuttl H98l [19921) . Ha can be observed in 
the o ptical in the local Universe (e.g., iGallego et al.1 
119951: ITresse et all 119981 : iKennicuttl 120081 ) . However, at 
z > 0.4, i t must be observed in the nea r infrared (e.g., 
Yan et all [1991 iGlazebrook et al.l Il999t IHopkins et atl 



spectroscopy (e.g., 



Hammer et al 1997; H ogg et alJl998t 



Gallego et atll2662tlTeplitz et al. 2003: R igopoulou et al 



2005J), or narrow-ba nd imaging (e.g 
2003t iLv et al.l I2007L hereafter Ly07; iTakahashi eTal 



Hippelein et al 



20001: ITresse et al.l l2002| ). or other less direct SFR indi- 



cators such as the [O II] AA3726,3729 doublet must be 
used. 

Compared to Ha, [O II] is only indirectly cou- 
pled to the ionizing continuum from massive stars, 
and is more sensitive to variations in metal abun- 
dance, excitation, and dust attenuation. Neverthe- 
less, because of its intrinsic strength and blue rest- 
frame wavelength, [O II] remains a good alternative 
SFR indicator for h i gh-redshift ga l axies (iKennicuttll 19981: 



Jansen et al.1 120011: iKewlev et al.1 12004 iMouhcine et al.1 
2005L iMoustakas fc Kennicuttll2006l ). In the past decade, 
a number of investigators have measured the star for- 
mation rate density, psfr, at high redshift by study- 
ing the [O II] luminosity function (LF), either using 
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20071 hereafter Takahashi07). Unfortunately, these stud- 
ies have been hampered by small sample size, small vol- 
ume probed, and an inconsistent treatment of dust ob- 
scuration. Methods like narrow-band imaging also suffer 
from difficulties in continuum subtraction and contami- 
nation from other emission lines. 

To circumvent these issues, we measure the [O II] lu- 
minosity function at z ~ 1 using data from the Deep 
Extragalactic Ev olutionary Probe 2 survey (DEEP2 2 ; 
iDavis et al1l2003D . DEEP2 has obtained high- resolution 
spectra for ~ 50, 000 objects over four separate fields, 
making it the largest existing spectroscopic redshift sur- 
vey of galaxies at these redshifts. We use these data to 
measure the [O II] luminosity function in four redshift 
bins at 0.75 < z < 1.45. 

In §2, we briefly describe our sample and the method 
used to calculate the LF. In §3 and §4 we present the 
observed [O II] LF and its evolution with redshift, re- 
spectively. Finally, in §5, we compute psfr in several 
redshift bins centered around z ~ 1, and we summarize 
our principal conclusions in §6. 

Throughout this work, we adopt a ACDM cosmology 
with fi m = 0.3, ft A = 0.7 and H = 70 km s" 1 Mpc" 1 . 
All magnitudes are on the AB system. 

2. DATA AND METHOD 

2.1. Data: DEEP2 DR3 

We select our sample from the DEEP2 third public 
data release (DR3), which includes BRI photometry and 
spectra for ~ 50, 000 galaxies in four widely separated 

2 http://deep.berkeley.edu 
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fields. The high-resolution (R = X/SX ~ 5000) spectra, 
which were acquired using the Keck-II/DEIMOS spectro- 
graph ijFaber et al.ll2003t) . span 6525 - 9120 A. We thus 
restrict our analysis to the redshift range 0.75 < z < 
1.45, where [O II] is measurable. To study the evolution 
of the [O II] LF, we further split each bin into four ad- 
ditional redshift bins: 0.75 < z < 0.93; 0.93 < z < 1.10; 
1.10 < z < 1.28; and 1 28 < z < 1.45 We refer 
the reader to IDavis et all (|2003f ) . ICoil et alJ (|2004f ). and 
IDavis et al.l (|2007ft for additional details regarding the 
DEEP2 survey. 

In addition to the flux cut, R < 24.1 mag, the DEEP2 
team applied the following color cuts to preselect galaxies 
at z > 0.7 in Fields 2, 3, and 4: 

B-i?<2.35(i?-/)-0.25 or; 
i? /> 1.15 or; 

B-R<0.50. (1) 

In the Extended Groth Strip (EGS, Field 1), to test their 
sample preselection method, the DEEP2 team did not 
apply the BRI color cuts. However, in order to make 
all four fields consistent, here we apply the color cuts 
given by equation ([T]) to the redshift and photometric 
catalogs in the EGS field. We apply additional angular 
cuts to avoid survey edges and to exclude gaps in the 
spectroscopy. Our final catalog contains 36, 118 objects, 
of which 24,729 h ave accurate redsh ifts (quality Q = 3 
or 4 as defined by IDavis et alJ l2007f ) . The areas of the 
four fields are 0.38, 0.56, 0.88 and 0.63 deg , respectively, 
and the total area is 2.45 deg 2 . 

Next we describe how we derive the integrated [O II] 
luminosity for each galaxy. Unfortunately, the DEEP2 
spectra are not flux-calibrated; therefore we infer the 
[O II] luminosity by multiplying the rest-frame emission- 
line equivalent width (EW) measured in the optical spec- 
tra, by the continuum luminosity around 3727 A inferred 
from fitting the BRI photometry. This method has the 
advantage that it is insensitive to the absolute calibra- 
tion of the optical spectra. However, the 1" wide slit in 
the DEEP 2 survey may not enclose all the line-emitting 
regions in the galaxy, thus this method does assume that 
the relative intensity of star formation and stellar light 
inside and outside the slit does not vary significantly. 
Without spatially resolved spectroscopy, however, this 
assumption is difficult to test directly. 

To measure the emission-line EW, we model each com- 
ponent of the [O II] AA3726, 3729 doublet simultaneously 
using two Gaussian profiles, constrained to have the same 
intrinsic velocity width and a fixed wavelength separa- 
tion, and use a smooth B-spline to estimate the contin- 
uum level around [O II] . Dividing the total [O II] flux by 
the detected continuum yields the rest- frame EW([0 II]) 
in A. We then fit the broad-band BRI photometry at the 
known redshift to obtain the best-fitting spectral energy 
distribution (SEP) using the deep_kcorre ct routine in 
kcorrect 3 (v4.1.4; iBianton fc Roweisll2007T ). Finally, we 
multiply the EW by the flux-density of the best-fitting 
SED at 3727 A to obtain the integrated [O II] luminos- 
ity. Because the effective wavelengths of all the BRI 
filter bandpasses lie blueward of 3700 A above z ~ 1.2, 
this technique does require extrapolating the best-fitting 

l! http : // cosmo .nyu. edu/blanton/kcorrect 



SED beyond the effective /-band wavelength to estimate 
the 3727 A continuum luminosity for the galaxies in our 
highest redshift bin; however, the uncertainty introduced 
by this extrapolation is negligible. 

Our [O II] measurement tech nique is simi lar to 

that used by the DEEP2 team (jWeiner et al.l 120071 : 
ICooper et al.l I2001I ). although the procedures used to 
compute if-corrections are totally independent. A com- 
parison of our measurements shows no systematic differ- 
ences and a ~ 25% scatter, which is comparable to the 
typical measurement error. We also checked our work by 
replacing our measurements with theirs, and obtained 
consistent results. In the following analysis we use our 
measurements . 

We consider an [O II] measurement with signal to noise 
ratio (S/N) larger than 5 as reliable. In detail, our re- 
sults are not sensitive to the specific S/N cut used since 
it only affects weak [O II] detections, for which we are 
incomplete anyway. For example, using S/N > 2 has no 
significant effect on our conclusions. Our final sample of 
[O II]-cmitting galaxies contains 13, 944 objects. 

2.2. Method: 1/Knax Method 

To calculate the luminosit y function, w e use the non- 
parametric 1/Vmax method (jFeltenl [T9761. For a given 
galaxy, we calculate: 

K iax = \j dnf™ dz d[D j z z)3] f(z), (2) 

which is suitable for a spatially flat universe. The angu- 
lar integral is limited to the DE EP2 a r ea, D r .{z) is the 
comoving distance as defined by Hogg (1999), and f(z) 
expresses the probability of selecting each galaxy in our 
sample. We assume that f(z) is given by the product 
of four quantities: / = /target x /success X / covor x / cut , 
where /target is the rate at which a source of a given R- 
band magnitude and B — R and R — I color was targeted 
[see eq. Q], /success is the rate at which a redshift was 
successfully measured, / CO ver is the fraction of DEEP2 
spectra where the wavelength coverage includes the red- 
shifted [O II] line, and / cu t is a step function, which wc 
determined using the BRI magnitudes and the [O II] 
S/N, as described below. We experimented with includ- 
ing the galaxy surface brightness as a fifth variab le in 
the completeness function (see, e.g.. lLin et al.ll2008l ). but 
found that it did not significantly affect our measured 
LFs. 

For a given gala xy, we calcula te Vmax using a Monte 
Carlo method fsee lBlantonll2006l ). We randomly choose 
1200 values of redshift z between z m - m and z max , uni- 
formly distributed in volume. For each redshift, we cal- 
culate what the magnitudes of the object would be us- 
ing deep_kcorrect. To estimate the [O II] S/N at each 
mock redshift, we first determine the mean noise spec- 
trum for every DEEP2 mask. Then, for a given galaxy, 
we randomly choose a mask and the corresponding noise 
spectrum and compare the noise at the actual observed 
wavelength of [O II] to that at the simulated wavelength 
as if the galaxy were at the simulated redshift, and obtain 
the new S/N. We then apply the appropriate flux cut, 
color cuts, [O II] S/N cut and the completeness function 
to determine the fraction of mock sources that would 
have passed our selection criteria. Finally, Vmax is given 
by the comoving volume multiplied by this fraction. 
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luminosity bin, we define the effective weight (W c s) by: 
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Fig. 1. — Observed [O II] luminosity function in four rodshift 
bins. The solid lines are the power-law fits to the bright part of the 
luminosity function. The dashed lines in the three highest redshift 
bins correspond to the solid line in the top-left panel. The three 
dashed-dotted lines correspond to three different faint-end slopes, 
—1.1, —1.3 and —1.5, and the dotted lines in each panel indicate 
the turnover in the luminosity function. 

When determining the completeness function, we as- 
sume a 100% redshift success rate for blue galaxies; 
that is, we assume that the only targeted blue galax- 
ies without successfully measured redshifts are those 
for which [O II] falls outside the wavelength range of 
the spectra, and thus outside the redshift range stud- 
ied here. He re, we define b l ue ga laxies using the same 
definition as IWillmer et al.l (|2006|) , who also made the 
same assumption regarding blue galaxies without well- 
measured redshifts. This assumption is especially rea- 
sonable in our analysis because a well-detected [O II] 
doublet will always result in a well-measured redshift. 
Finally, we tested our completeness function by calcu- 
latin g the .B-band LF for all galaxies and comparing 
with IWillmer etail (|2006l) . and found excellent statis- 
tical agreement, and no systematic differences. 

3. THE OBSERVED [O II] LUMINOSITY FUNCTION 

We measure the [O II] luminosity function in four 
redshift bins: 0.75 < z < 0.93; 0.93 < z < 1.10; 
1.10 < z < 1.28; and 1.28 < z < 1.45, according to 
where the [O II] doublet could be measured reliably in 
the DEEP2 spectra. The median redshifts of the [O II]- 
emitting galaxies within these four redshift bins are 0.84, 
1.01, 1.19, and 1.35, containing 5471, 3771, 3018 and 
1684 galaxies in each bin, respectively. 

The resulting [O II] LFs are tabulated in Table [1] and 
illustrated in Figure [TJ The error bars shown in the figure 
are 84.13% confidence Poisson upper limits and lower 
li mits estimated using the approximate formulas given 
bv lGehrelsl (fl986l) [ see eqs. (10) and (14) in that paper]. 
The distribution of the LF must follow a scaled Poisson 
distribution. To determine the scaling factor, within each 



E 



(Knax)j 



and the effective number (N e g 

1 



N r 



E 




(3) 



(4) 



We then calculate the upper limit and lower limit for 
the effective number, and multiply them by the effec- 
tive weight to obtain the upper limit and lower limit for 
the luminosity density in each bin. For comparison, the 
square root Poisson error gives: 
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which is commonly used in the literature. The Poisson 
errors, however, do not include the effects of cosmic vari- 
ance. Because we have four widely separated fields, we 
determine the error due to cosmic variance, a cv , by cal- 
culating the variance among the four independent fields, 

and list the results in Table [JJ 

Exa mining Figure (TJ the commonly used iSchechter] 
(|1976|) function is clearly a poor representation of the 
data. Instead, we model the observed [O II] LF in each 
redshift bin as a power law: 

$(logL)d(logL) = 10 (Q+1)(logi " 42 - 5)+/3 d(logL), (6) 

where L is £[oii] in erg s^ 1 , and a and /3 are dimen- 
sionless parameters. We find the best fitting parameters 
(a, 0) using a non-linear least square fit to the [O II] LFs 
weighted by the average Poisson errors. These parame- 
ters are presented in Table [2 Our results show that the 
bright part of each LF can be represented by a power law 
dN/dL oc L a with slope a ~ —3. The slope for the high- 
est redshift bin is the flattest: —2.79 ± 0.12. However, it 
is possible that the slope in this bin may be underesti- 
mated due to incompleteness (see below). 

Unfortunately, we are unable to constrain the faint 
end of the [O II] LF, especially in the highest redshift 
bin. Because we do not know the intrinsic number of 
[O II] emitters within a certain [O II] luminosity bin, we 
can not calculate the completeness of [O II] luminosity 
directly, but can only infer that from the completeness 
of broad-band properties. To analyze the completeness, 
and in particular to test the significance of the observed 
turnover (TO) in the LF (vertical dotted lines; Fig. [JJ), 
we construct two diagnostic diagrams. In Figure [2] we 
show the U — B versus Mb color-magnitude diagram 
for our sample. We plot the approximate i?-band com- 
pleteness limits for blue galaxies as dash ed lines, where 
the slo ping color cut has been defined by IWillmer et al.l 
(|2006l ). and the thick vertical dashed line roughly cor- 
responds to where the data and the color cut begin to 
deviate. In Figure [3] we plot the density distribution of 
points in the M B — log(L[ n]) plane in grey scale; the 
two contours enclose 50% and 80% of the points, respec- 
tively. The thin dashed lines in this figure all have a 
slope of —3.1 mag dex -1 . This slope is formally consis- 
tent with performing a mean ordinary least square fit to 
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Fig. 2. — Rest-frame color-magnitude diagrams. The thin dashed 
lines show the divis ion between blue galaxies and red galaxies (see 
IWillmc r ct al. 20QH), and the thick dashed lines show the approxi- 
mate Mb completeness limit for our sample. 

the data (|Isobe et al.lll990f ). For our purposes, we note 
that the line roughly bisects the distribution of points in 
each panel, i.e., it approximates the median relation be- 
tween i[on] and Mb- The thick horizontal dashed lines 
in each panel are equivalent to the thick vertical dashed 
line plotted in the respective panels in Figure El and the 
vertical dotted lines give the position of the turnover in 
the respective [O II] LF (Fig. [T]). The majority of the 
galaxies missing from our sample should be along and to 
the left of the thin dashed line. 

Figure [3] demonstrates that, brighter than the 
turnover, we expect the sample to be complete because 
the majority of the unobserved galaxies below the Mb 
completeness limit are along and to the left of the thin 
dashed line. Fainter than the turnover, the sample be- 
comes increasingly incomplete. In the two lowest red- 
shift bins, the turnovers appear to be significant, because 
fainter than the turnover the difference between the mea- 
sured LF and the extrapolation of the power law fitted 
to the bright part of the LF is so significant that it is 
unlikely to be due to incompleteness. In the two high- 
est redshift bins, however, because the Mb completeness 
limit is very bright, it is possible that the turnovers are 
artificial, caused by the incompleteness of the survey. In 
the highest redshift bin, the Mb completeness limit is so 
bright that the slope of the LF may be even steeper than 
we have derived. 

To summarize, we are unable to constrain the faint 
end of the [O II] luminosity function. However, we em- 
phasize that the evolutionary analysis presented in §4 
is unaffected by our inability to constrain the faint-end 
slope, because we restrict our analysis to the bright part 
of the [O II] LF where we are statistically complete. In 
§5, when integrating the [O II] LF to obtain an estimate 
of the SFR density, we do make some simplified assump- 
tions regarding the form of the faint end of the LF. 



Fig. 3. — Mb — l°g(^[on]) distribution. The dotted lines are 
the same as in Figure [T] The thick dashed lines show the Mb 
approximate completeness limit for blue galaxies (see Figure [2{ . 
The thin dashed lines have a slope of —3.1 mag dex~ x , and show 
the approximate median relation between Mb and log(L[on])- 

In Figure [¥J we compare our results at z ~ 1.2 
with three other recent coeval measurements of the 
[O II] LF based on narrow-band observations (Ly07, 
Takahashi07). Using Subaru/NB816 narrow-band imag- 
ing of the 875 arcmin 2 Subaru Deep Field (SDF) by 
iKashikawa et all (|2004r i, Ly07 and Takahashi07 identi- 
fied 894 and 602 [O II]-emitting galaxies, respectively. 
In the 2 deg 2 Cosm ic Evolution Survey field (COSMOS; 
IScoville et atll2007| ). Takahashi07 used the same narrow- 
band filter to identify 5824 [O Hi-emitting galaxies at 
z ~ 1.2 (jTaniguchi et all l2007| ). Over the range of 
[O II] luminosity where all the surveys are complete, 
10 42 erg s" 1 < i^on] < 10 42 7 erg s" 1 , we find that 
our LF from DEEP2 agrees very well with the luminos- 
ity functions derived by Ly07 and Takahashi07 in the 
SDF and COSMOS fields, respectively. It is not clear 
why the LF of the SDF derived by Takahashi07 is so dis- 
crepant with the other surveys, although cosmic variance 
may play a role. Nevertheless, this comparison shows 
that: (1) our assumptions regarding the faint-end slope 
are reasonable; and (2) the bright end of the [O II] LF is 
clearly a power law, not a Schechter function. 

4. THE EVOLUTION OF THE [O II] LUMINOSITY 
FUNCTION 

The uncertainty in the faint end of the [O II] LF pre- 
vents us from obtaining a reliable estimate of the total 
[O II] luminosity density. Therefore, we focus instead 
on the integrated number density, <f> = J L Q(L) dL, of 
the strongest [O II]-emitting galaxies. Assuming that 
the turnover luminosity in the highest redshift bin is 
< 10 42 erg s _1 , we integrate our model of the LF in 
each redshift bin over I/foii] > 10 42 erg s _1 , and present 
the results in the top panel of Figure [5] and in Table [5J 
The horizontal error bars indicate the range of each red- 
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Fig. 4. — Comparison of our [O II] luminosity function at 
2 ~ 1 . 20 with other recent measurements from the literature based 
on narrow-band imaging. The open triangles and the solid line rep- 
resent our data and best-fitting power law (Fig.[T] lower-left panel). 
The open squ ares represent th e LF in the Subaru Deep Field (SDF) 
calculated bv lLv et all 12007T ) atz~ 1.18. The open and filled cir- 
cles correspond to observations of galaxies at z ~ 1.19 in the SDF 
and COSMOS fields, respectively lITakahashi et al.|[200^ . Finally, 
the shaded region represents our assumption for the faint-end slope 
when we convert the [O II] LF to SFR density in §5. 

shift bin, and the vertical error bars are given by the 
cosmic variance among the four fields, which dominate 
the error budget. We find that the total number density 
of the strongest [O II]-emitting galaxies, i.e., those with 
iron] > 10 42 erg s _1 , declines by a factor of > 2.5 be- 
tween z ~ 1.35 and z ~ 0.84. A linear fit to the four 
points gives a slope of ~ 0.9 dex per unit redshift: 



log [<t>(L [01I] > 10 42 erg s" 



a z 



b, 



(7) 



with (a, 6) = (0.90 ±0.14, -4.43 ±0.18) and (f>(L [olI] > 
10 42 erg s _1 ) in Mpc -3 . This fit is shown as the dashed 
line in Figure [5j 

The observed decrease in the number density of [O II]- 
emitting galaxies may be caused by the decline of the 
overall number density of galaxies in the Universe, or 
by the evolution of the luminosity function itself. We 
attempt to constrain the amount of evolution by mea- 
surin g the luminosity a t a fixed space density (see, 
e.g., iBrown et al.l 120071 ). We define a characteris- 
tic luminosity, ijon]' wnere the space density equals 
-LQ-3.5 dcx _1 Mpc -3 , which has the advantage that it is 
independent of the faint end of the LF. We show £f OII i 
as a function of redshift in the bottom panel of Figure 03 
and list the results in Table [2j Once again, the vertical 
error bars arc dominated by cosmic variance. We find 
that ipii] declines by a factor of ~ 1.8 between z ~ 1.35 
and z ~ 0.84. We perform a linear fit to the four points 
and obtain a slope of ~ 0.5 dex per unit redshift: 



logL 



[Oil] 



c z + d, 



(8) 



Fig. 5. — Evolution of the [O II] luminosity function. (Top) Evo- 
lution of the total number density in Mpc -3 of the strongest [O II]- 
cmitting galaxies with I/[on] > 10 42 erg s — 1 . The dashed line is 
the linear fit to the data given by eq. (7) . (Bottom) Evolution of the 
characteristic luminosity, in erg s — 1 , defined by where the space 
density of [O II]-emitting galaxies equals 10 -3 ' 5 dex - iMpc - 3 . 
The dotted line is the linear fit to the data given by eq. (8). 

with (c,d) = (0.53 ± 0.07,41.79 ± 0.09) and L [0 n] in 
erg s _1 . The resulting fit is shown as the dotted line 
in Figure [5] 

To summarize, we find that the total number density 
of the strongest [O II] -emitting galaxies has declined by 
a factor of > 2.5 between z ~ 1.35 and z ~ 0.84, when 
the Universe aged from 4.6 Gyr to 6.4 Gyr. This de- 
cline may be driven by a decline in the overall number 
density of galaxies in the Universe, or by a fading of 
the [O II] LF. Unfortunately, we are unable to establish 
whether number-density evolution, luminosity evolution, 
or a combination of both is responsible for the observed 
evolution. Nevertheless, if we assume that the observed 
change in the [O II] LF is predominantly due to luminos- 
ity evolution, that L[oii] is proportional to the SFR (see 
§5), and that the typical dust attenuation in galaxies at 
fixed [O II] luminosity does not change significantly over 
0.75 < z < 1.45, then this result implies that the SFR in 
galaxies declines by ~ 25% per Gyr during this epoch. 
In the next section we integrate the full [O II] LF using 
some simple assumptions, and compare our results with 
other estimates of the SFR density at z ~ 1. 

5. THE EVOLUTION OF THE STAR FORMATION RATE 
DENSITY 

Converting the observed [O II] luminosity into a SFR 
is subject to considerable random and systematic un- 
certainties, arising from variations in dust attenuation, 
metallicity, and excitation among s t ax-forming galaxies 
(lKennicuttJll992tlJansen et alJl2001tlKewlev et aij |2004t 
iMouhcine et al.ll2005t iMoustakas et al.ll2006[) . Neverthe- 
less, our measurement of the [O II] luminosity function 
affords a valuable opportunity to constrain the SFR den- 
sity of the Universe during an important epoch of cosmic 
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Fig. 6. — SFR density, PsfRj v ersus redshift based on various 
multi-wavelength SFR indicators (Hopkins 2004). Our four esti- 
mates of psFR are shown as large filled squares, where the vertical 
error bars are obtained by allowing the turnover positions and the 
slope of the faint end of the [O II] luminosity function to vary in a 
sensible way (see text for details) . 

history. 

In a recent analysis, iMoustakas et al.l (|2006T > showed 
that dust reddening, as derived using the Ha/H/3 Balmcr 
decrement, is responsible for the bulk of the scatter in 
[O II] as a SFR indicator, while variations in metallicity 
and excitation are second-order effects for most galaxies. 
Unfortunately, the DEEP2 spectra do not span a suffi- 
ciently wide wavelength range to include Ha, H/3, and 
other emission-line diagnostics of the metallicity and ex- 
citation. Therefore, we use the e mpirical correlation de- 
rived by IMoustakas et al.l (|2006[ ) between the absolute 
.B-band magnitude, and the L[ n]/SFR ratio. This cal- 
ibration statistically accounts for the gross systematic 
effects of reddening, metallicity, and excitation, all of 
which correlate with optical luminosity, and has been 
shown to work rea sonably well for star-fo r ming galaxies 
at 0. 7 < z < 1.4 (jMoustakas et al.l 120061: ICooper et all 
12001 . Note that an [O II] SFR conversion that is in- 
dependent of luminosity (e.g., iKennicuttl 119981) would 
severely underestimate the SFR, because luminous, star- 
forming galaxies tend t o be dustier and more metal-rich 
(jMoustakas et al.ll2006l) . 

Another potential concern is that some fraction of the 
[O II] emission might be arising from an active galactic 
nucleus (AGN) rather than star fo rmation. Although the 
AGN - sensitive [N Hi /Ha ratio (jVeilleux fc Osterbro"ckl 
Il987t iKewlev et all 1200 it ) lies in the near-infrared at 
z > 0.75, AGN that contribute significantly to the op- 
tical emission-line spectrum also can be identified using 
[O III] A5007/H/3. In the DEEP2 spectra, [O III] /H/3 is 
measurable for galaxies at 0.75 < z < 0.85, comprising 
roughly one-third of the sample in our lowest redshift 
bin. Among these objects, we find 2% (32/1485) with 



log([0 III]/H/3) > 0.7, indicative of AGN activity. Wc 
can also leverage deep X-ray obser vations of the EGS 
to identify AGN (jLaird et al.ll2009t) . Among the 2579 
galaxies in this field, we find that 1% (28/2579) are also 
X-ray point sources. These results reveal that power- 
ful AGN constitute a negligible fraction of the sources 
in our sample. However, even if we have significantly 
underestimated the fraction of AGN in our sample, de- 
tailed studies show that the physical conditions in the 
narrow-line regions of AGN in the local Universe dis- 
favor [O II] emission, which is o ne advantage of using 
[O II] as a SFR tracer ()Holl2005l ). Indeed, we will show 
below that our [O II]-based estimate of the SFR density 
at 0.84 < z < 1.35 agrees remarkably well with other 
multi-wavelength studies. Therefore, we conclude that 
AGN contamination is a negligible source of error on our 
results. 

Before integrating the observed [O II] LF to derive the 
SFR density, psfr, we must make some assumptions re- 
garding the form of the faint end of the LF (see §3). 
First, we allow the luminosity of the turnover in the LF 
in each redshift bin to vary over a sensible range of val- 
ues to account for the uncertainties in our completeness. 
Specifically, for the two lowest redshift bins, we assume 
log(£ T o) = 41.68 ± 0.10 and 41.75 ± 0.10 erg s _1 , re- 
spectively, while for the two highest redshift bins, wc 
adopt a fainter lower limit: log(I/To) = 41.93^20 an d 
42.001°;^ erg s" 1 . 

Second, we must assume a form for the [O II] LF fainter 
than the turnover luminosity. P revious studies (e.g. 
iGallego et all l2002t iLv et al.l l2007f ) have assumed that 
the [O II] LF is a Schechter function, which is a power 
law at the faint end. However, the faint-end slope, af a i n t, 
is usually not well-constrained. Consequently, hereafter 
we allow af a i n t to vary between —1.1 and —1.5, which 
brackets the value, af a i n t = — 1.3 ± 0.2, that we measure 
from the lowes t redsh ift bin in Figure [TJ For comparison, 
IWillmer et al.l (|2006f ) assumed a faint-end slope of —1.3 
for the i3-band LF of blue galaxies at 0.3 < z < 1.3. 

To summarize, we assume that the [O II] luminosity 
function is a double power law with slope: 

f af a i n t = -1.3 ± 0.2 for L [ou] < L TO 
\ a for L [OII ] > L TO 

where a is derived from our fit to the bright part of the 
LF where we are complete. Given these assumptions, we 
integrate the observed LFs and list the results in Tabled 

Instead of using the S-band luminosity of each individ- 
ual galaxy, wc obtain a statistical estimate of Mb for each 
object from £[oii] using the thin dashed line in Figure EH 
We then calculate the appropriat e L[pn|/SFR conv ersion 
factor by interpolating Table 2 in Mousta kas et al.1 (|2006l ) 
to derive the SFR. 

Figure [5] compares our measurements of psfr at z ~ 1 
against a large compilation of multi- wavelen gth measure- 
ments from the literature bv lHopkinsl (|2004h . Our results 
agree remarkably well with these independent measure- 
ments considering the uncertainties in converting 
into a SFR, and our incompleteness at the faint end of 
the LF. 

6. CONCLUSIONS 

Because its blue rest-frame wavelength and intrinsic 
strength allow it to be measured up to z ~ 1.6 in the op- 
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tical, the [O II] AA3726, 3729 doublet plays a unique role 
in the study of galaxy evolution. We have used spec- 
troscopy of ~ 14, 000 galaxies from the DEEP2 galaxy 
rcdshift survey to measure the [O II] luminosity function 
at 0.75 < z < 1.45. Our sample is orders-of-magnitude 
larger than previous spectroscopic studies, over a con- 
siderable larger area spanning four independent fields, 
allowing us to minimize the systematic effects of cosmic 
variance. Our principal results are given in Tables [1] and 
[21 and illustrated in Figures Q] and [5] We found that the 
bright part of the [O II] LF is well-represented by a power 
law dN/dL oc L a with slope a ~ —3. However, survey 
incompleteness prevented us from constraining the faint 
end of the LF. 

We measured the evolution of the [O II] LF using 
two quantities that only rely on the bright part of the 
[O II] LF where we are statistically complete. First, 
we calculated the total number density of galaxies with 
L[ou] > 10 42 erg s -1 , and found that it has declined by 
a factor of > 2.5 between z ~ 1.35 and z ~ 0.84. Second, 
we calculated the characteristic luminosity, the luminos- 
ity where the space density of [O II]-emitting galaxies 
equals 10 -3 ' 5 dex _1 Mpc -3 , and found that it has de- 
clined by a factor of ~ 1.8 over the same redshift inter- 
val. Assuming that the [O II] luminosity is proportional 
to the SFR, these results imply that the SFR in galaxies 
declined by ~ 25% per Gyr during this epoch. 

Finally, we used the empir ical calibration b etwee n 
L[on] and SFR published by iMoustakas et al.l (|2006() . 
and adopted some simple assumptions regarding the faint 
end of the [O II] LF, to obtain an estimate of the inte- 
grated SFR density, psfr, in four redshift bins centered 
around z ~ 1. We found that, despite the considerable 
uncertainties, the evolution we measure is consistent with 



previous measurements based on a variety of indepen- 
dent, multi-wavelength SFR indicators. 
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TABLE 1 

OBSERVED [O II] LUMINOSITY FUNCTION 



log £[OII] 
(ergs s ) 


0.75 < z < 0.93 


0.93 < z < 1.10 


1.10 < z < 1.28 


1.28 < z < 1.45 


* CTcv N gal 


* CTcv AT Ba! 




* O-cv N ga i 



41.15 
41.30 
41.45 
41.60 
41.75 
41.90 
42.05 
42.20 
42.35 
42.50 
42.65 
42.80 
42.95 
43.10 
43.25 
43.40 
43.55 
43.70 
43.85 



49.49 1 



46.93_ 187 

-1.37 
-1.31 

l'-t>8_0.90 
,+0.80 
-0.72 
-, + 0.55 
-0.50 
-,+0.35 
-0.29 
-,+0.49 
-0.35 
2+0.15 
-0.09 
2+0.21 
-0.11 
7+0.15 
-0.10 
□+0.12 
-0.06 



29. 00 1 



7.91] 
4.70^ 
1.62j 
l.lOj 
0.18j 
0.18j 
0.27j 
0.09^ 



0.04 



+0.08 



7.30 



553 



+3.58 
-3.34 

60.57±3 2 7 19.41 918 

56 - 07 -2.ii 15 - 88 1095 

194 -7 RA 101 g 



7.64 

5.82 669 

3.68 429 

1.57 176 

1.63 104 

0.55 35 

0.44 19 

0.11 4 

0.24 3 

0.12 7 

0.09 2 

0.08 1 



20.69+J ; 47 



28.25 
35.57 
31.30 
20.38 
11.60 
6.14 
2.77 
1.13 
0.78 
0.27 
0.07 
0.15 
0.06 



+2.03 
1.90 
1.81 
1.72 
1.44 
1.38 
+ 0.93 
0.89 
+0.72 
0.68 
+0.81 
0.72 
+0.60 
0.50 
+ 0.54 
0.38 
+0.22 
0.17 
+0.14 
0.10 
+0.09 
0.04 
+0.12 
0.07 
+0.09 
-0.04 



0.03 
0.03 



+0.08 
0.03 



+0.08 
0.03 
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.08 
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1 



18.8611$ 
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-,+0.93 
-0.88 
q co+0.60 

,+0.42 
-0.39 
,+0.33 
-0.29 
-+0.24 
-0.20 
7+0.12 
-0.09 
4+0.18 
-0.12 
7+0.11 
-0.07 

jj+0.08 

-0.04 
,+0.08 
-0.03 
,+0.06 
-0.02 



16.60"* 



4.52] 
2.33^ 
l.lSj 
0.27j 
0.34j 
0.17j 
0.06j 
0.03j 
0.03^ 



1. 


57 


349 


2 


03 
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1. 
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658 





11 


506 





96 
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33 


37 





11 


10 





21 


9 





10 
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06 


2 





05 


1 





04 


1 



12.17 



-1.19 
-1.09 
-1.18 



2.22 207 
1.74 303 
376 



1 5 5 £ }~ t ~ 

16.66l};j5 °- 94 

9.70l°; 7 i 0.53 247 

5.59lo;4 7 0.97 155 

3.99I045 0.38 108 

n + 0.'32 

-0.28 
7+0.23 



1.93] 



0.48 



0.87Iois °- 28 
55+° :i9 

u - oo -0.14 

o n+ 011 

u - ii -0.06 
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1-0.08 
-0.03 



0.03] 



52 
24 

0.18 15 
0.12 3 
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0.08 4 
0.24 3 
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Note. — $ is in units of 10 dex Mpc , a cv is the la uncertainty in 3> due to cosmic variance, and N ga i is the number of galaxies in 
each bin. For reference, the median rcdshifts of the sources in each of the four redshift bins arc 0.837, 1.005, 1.191, and 1.349, respectively. 
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TABLE 2 

PARAMETERS OF THE OBSERVED [O II] LUMINOSITY FUNCTIONS 



Quantity 


0.75 < z < 0.93 


0.93 < z < 1.10 


1.10 < 2 < 1.28 


1.28 < z < 1.45 


^■median 


0.837 


1.005 


1.191 


1.349 




5471 


3771 


3018 


1684 


Age of Universe (Gyr) 


6.46 


5.73 


5.06 


4.59 


a 


-3.01 ±0.07 


-2.98 ± 0.08 


-3.03 ± 0.07 


-2.79 ±0.12 





-4.01 ±0.04 


-3.86 ±0.04 


-3.67 ±0.03 


-3.49 ±0.03 


log*(L [0 ii] > 10 42 ergs- 1 ) 3 


_o fi7 +o.io 


, c-o+0.07 
°- oo -0.08 


n oq + 0.08 

°' oz -0.10 


o 91 +0.01 
°" il -0.01 


logL T o c 


42.24 ± 0.06 


42.32 ± 0.04 


42.42 ± 0.03 


42.50 ±0.01 


41.68 ±0.10 


41.75 ±0.10 


41 QQ+0.10 


42.00_°; ° 


log p(L[Oii]) d 


39.34_°;» 


39.40±g:ig 


39.44+°?° 


39.47_°; 24 


PSFR C 


0.091°;°?, 


0.10_°;°I 


12+°'° 7 
u - lz -0.03 


o.i4±°;°I 



a <j> = $(L) dL is in MpcT 3 . 

b The characteristic luminosity, LT on ,, where the space density equals 10 _3 ' 5 dex _1 Mpc -3 , in 



erg s 

c Luminosity of the turnover in the luminosity function, Xto, m er g s 
d Integrated luminosity density, p(L[on]), in erg s" 1 Mpc~ 3 . 
6 Star formation rate density in M@ yr _1 Mpc -3 . 



